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ABSTRACT 

The abundance stratification in the ejecta of the normal Type la Supernova 2002bo 
is derived fitting a series of spectra obtained at close time intervals. A Montecarlo code, 
modified to include abundance stratification, is used to compute synthetic spectra at 13 
epochs in the photosphcric phase, starting 13 days before B maximum. A description 
of the abundance distribution above 7600 km s~^ is thus obtained. Abundances in 
deeper layers, down to zero velocity, are derived from models of two nebular-phase 
spectra. Elements synthesised in different stages of burning are significantly, but not 
completely mixed in the ejecta. A total ^^Ni mass of 0.52 M0 is derived. Evidence for 
intermediate-mass elements at high velocities (> 18, OOOkms"^) is found, most clearly 
in Si II 6355 A, Can H&K, and in the Call IR triplet. Carbon lines are not seen at any 
velocity, with possible implications on the progenitor/explosion scenario. A synthetic 
bolometric light curve computed using the inferred abundance distribution is in very 
good agreement with the observed one, providing an independent check. In particular, 
the fast rise of the light curve is reproduced very well. This is due to outward mixing 
of ^^Ni, which is clearly well determined by the spectroscopic modelling. 
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1 INTRODUCTION 



It is widely agreed that Type la Supernovae (SNe la) result 
from the thermonuclear explosion of Carbon-Oxygen white 
dwarfs (WD) in binary systems. Two possible progenitor 
configurations are usually considered. In the single degener- 
ate scenario a WD accretes mass from a Roche lobe-filling 
Red Giant companion. When the WD reaches a mass close 
to the Chandrasekhar limit (« 1.38 Mq), carbon burning is 
triggered by compressional heating n ear the centre. After a 
few thousand years of quiet burning ilben fc T\itu kovlll984t 
IWebbinkllTosi IWooslev. Wunsch fc Kuhlenll2004^ . a ther- 
monuclear runaway occurs which disrupts the star. In the 
double degenerate scenario, two low-mass WDs in a close bi- 
nary system with total mass exceeding the Chandrasekhar 
mass lose angular momentum via the emission of gravita- 
tional waves, which ultimately leads to the mergin g of the 
two stars and to a thermonuclear runaway (Whel an fc IbenI 
I1973I: lNomotolll982l: iHan fc Podsiadlowskill2004l) . Potential 
progenitor systems for both channels have been detected, 
but their numbers ar e too low to explain the frequency of 
occurrence of SNe la iCappellaro. Evans fc Turattdll999^ . 



Although the apparent predictability and homogene- 
ity of SNe la, together with their brightness, has mo- 
tivated their use as standardisable co smological candles 
jPerlmutter et alJl999l : lRiess et alll998l) . much remains un- 
known about their properties, as regards both global physi- 
cal processes and the peculiarities of individual objects. 

For example, the details of the explosion process are 
still unclear. Ignition is supposed to start near the centre. A 
subsonic (deflagration) wave, often called a "flame", trav- 
els outwards, burning part of the WD to nuclear statis- 
tical equilibrium (NSE). The subsonic speed of the defla- 
gration front prevents the WD from being burned to NSE 
entirely. Partial burning of C and O results in the produc- 
tion of intermediate mass elements (IME): Si, S, Mg and 
Ca dominate the ear ly spectra of SNe la. The prompt det- 
onation mechanism llArnettiligei) . on the other hand, is 
inconsistent with the spectra, as it fails to produce suffl- 
cient amounts of IME. There is still no agreement as to 
whether the explosion continues as a deflagration, becom- 
ing strongly turbulent (Nomoto, Thielemann & Yok oi 1984; 
I Wooslev. Axelrod fc Weavejll984l: iNiemever & Hille brandtl 
I1995I: iReinecke. Hillebrandt fc NiemeveJ 120021) or it turns 
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to a supersonic detonation. Recent 3D deflagration models 
leave significant amounts of unburned mater ial both in the 
outer an d in the inner parts o f the electa llGamezo et alJ 
12003 : Reinecke et al . 1999: Tr avaglio et alJl2004) . and pro- 
duce too little radioactive ^^Ni. This could affect the light 
curve of SNe la, which is powered by ^^Ni through the en- 
ergy released in the decay chain to s table ^^Fe via ^^Co. 
Presentl y, delayed detonation models jHoflich fc Khokhlovl 
[l996; Iw amoto et al. 1999') describe better the abundances 
of Fc-group elements and IME at high velocities, as in the 
detonation phase burning is boosted in the outer regions of 
the envelope. These models also have the flexibility to pro- 
duce different amounts of ^^Ni and IME, but this depends 
on the ad-hoc assumption of the time and position of occur- 
rence of the Deflagration to Detonation Transition (DDT, 
lKhokhlov|[l99ltlWooslev fc Weaveilll994) . 

Another area that needs clarification is the final dis- 
tribution of the elements in the ejecta. Light curve stud- 
ies typically make use of the direct results of the explosion 
models, where the abundances are stratified, but spectral 
analysis has favoured mixing in the ejecta, at least above 
some velocity l^ranch. et al. 19^rIi . Mixing was fav oured 
over stratification in an LTE study ljHarknesslll'99lh . but 
no further investigations have been performed, although the 
results could be very useful to discriminate between different 
scenarios for the explosion. 

The European Research and Training Network (RTN) 
"The Physics of Type la Supernova Explosions" was set up 
to study SNe la by means of very good time series of multi- 
wavelength observations and detailed models of a sample of 
nearby objects. SN 2002bo was its first target iBenetti et alJ 
|2004) . In this paper we model the spectra and light curve of 
SN 2002bo, and attempt to extract information about the 
abundance stratification in the SN ejecta. 

In Section |21 we summarise the basic properties of 
SN 2002bo. The method of the analysis is discussed in Sec- 
tion Models of the photospheric epoch spectra are pre- 
sented in Section|l] while Section|5]focuses on models of the 
nebular spectra. The derived abundance distribution is pre- 
sented and discussed in Sectional In Section a synthetic 
light curve computed on the basis of the derived abundance 
distribution is presented and discussed. In Section |H| we dis- 
cuss the behaviour of the line velocities, and in Section 9 the 
results are summarised. 



2 OBSERVATIONS 

SN 200 2bo was disco vered independently by Cacella and 
Hirose llCacella et alJ 12002) on UT March 9.08 and 9.505, 
2002. It is located in NGC 3190 at a distance modulus 
U = 31.67 fecn otti ct al. 2004, Ho = 65kms~^ Mpc'^). 
llBenetti et alTll2[)04ll) present photometry and spectroscopy 
covering from ~ 13 days before to 368 days after the esti- 
mated epoch of B maximum. They quote a B-band risetime 
of 17.9 ± 0.5 d, a decline rate Ami5(-B) = 1.13 ± 0.05, 
a reddening E{B — V)obs = 0.43 ± 0.10 and a reddening- 
correct ed Mb = -19.41 ± 0.42, making SN 2002bo a normal 
SN la iBranch. Fisher Nugendligga^ . However, compari- 
son with other SNe la reveals peculiarities, e.g. with respect 
to line velocities. 



3 METHOD OF ANALYSIS 

In order to describe the early-time evolution of a SN spec- 
trum (and in particular of a SN la) it is in principle neces- 
sary to solve the time-dependent problem of 7-ray transport, 
heating of the ejecta gas, and transport of the radiation. 
Ideally, this should be done in NLTE, and possibly in 3D to 
account for deviations from spherical symmetry which are 
probably not rare occurrences in SNe la. 

This is a complex problem, and it heavily relies on using 
a predefined model of the explosion. While the development 
of a code as briefly outlined above is under way, to meet 
also the requirements of the 3D explosion calculations which 
are becoming available from various groups, we nevertheless 
want to obtain as much information as possible from the 
observational data sets that are already available. 

SNe have the unique feature that they unfold their con- 
tent before our eyes: as they expand, deeper and deeper 
layers are exposed and contribute to line formation. Fitting 
a closely knit time series of spectra offers therefore the op- 
portunity to study the properties of the ejecta as a function 
of depth, in what may be viewed as a sort of CAT scan. Very 
early observations, covering the first week after the explo- 
sion, are particularly important. At these epochs the ejecta 
at velocities > 18,000kms~^ are dense enough to produce 
spectral lines. Therefore it is possible to derive information 
about the material at the surface of the WD. Also, it may 
be possible to detect signatures of an interaction between 
the ejecta and circumstellar material, which may originate 
from the WD or the companion star, thus shedding light on 
the properties of the progenitor system. The detail and the 
amount of information that can be derived from this pro- 
cedure clearly depends on the number of spectra that are 
available, and on the size of the time steps. Physical condi- 
tions in the first 8-10 days after the explosion change very 
rapidly compared to later epochs. Also, information about 
both the progenitor and the possible interaction with cir- 
cumstellar environment is present in the outermost layers of 
the ejecta. Therefore, it is desirable to have frequent very 
early observations, while a lower frequency is sufficient near 
maximum light and later. SN 2002bo has early and frequent 
data, and is therefore very well suited for this type of study. 

The best way to model the spectra independently of 
a particular explosion model is to use a simple but flexi- 
ble code, that can be adapted to extract information from 
the spectra: we base model calculations in the photospheric 
epoch o n a Montecarlo (MC) code developed in a se r ies of 
pape r s llAbbott fc Lucvlll985l: iMazzaH fc Lucvlll993l : ILucvI 
I1999I: 'Mazzali' '2000'), and successfully applied to various 

SNe la (Mazzali et al 1993. : .Mazzali. Danziger fc Turattol 

ll995l:lMazzaH et al boOSal) . 

The code uses a sharp lower boundary, below which all 
energy is assumed to be deposited, and follows the propaga- 
tion of energy packets in a spherically symmetric envelope. 
The density in the envelope depends on a model of the ex- 
plosion, but abundances can be arbitrarily chosen to fit the 
spectra. lonisation and excitation conditions are computed 
using a modified nebular approximation, w hich was found to 
be a good approximation to NLTE results iPauldrach et all 
Il997fl . Energy packets can undergo electron scattering and 
line absorption followed by reemission. The latter is treated 
using the Sobolev approximation, which is appropriate for a 
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SN envelope, and the process of photon branching. Finally, 
the emerg ent spectru m is computed using a formal integral 
approach l)Lucj l999'). 

For the calculations in this paper, we have adopt ed the 
W7 density- velocity structure jNomoto et al.lll984) . This 
profile is similar in most computed models, whether defla- 
grations or delayed detonations. Other inputs required by 
the code are the photospheric velocity Vph, the emergent lu- 
minosity L, the epoch t^xp (time since explosion), and the 
abundances. In order to study the abundance distribution in 
the ejecta, we have modified the code to allow for radially 
varying abundances. The detailed procedure is as follows. 

An abundance structure is set up where the shells are 
the photospheric radii of the observed spectra. Therefore, 
the position of this additional grid is not fixed a priori but 
is set during the model calculations. These start with the 
earliest spectrum, since it is the one with the highest pho- 
tospheric velocity. A homogenous abundance distribution 
above this first photosphere may be assumed, or a number 
of supplementary shells (typically ^ 3) may be introduced 
at higher velocities. An optimal model is then computed 
and its chemical composition, as well as the velocities of the 
boundaries of the abundance shells, are stored. At the epoch 
of the next spectrum in the time series the photospheric 
velocity has decreased, and deeper parts of the ejecta con- 
tribute to forming the spectrum. Although some conditions 
(density, temperature, etc.) have changed in our previously 
considered shells, the relative abundances have not, except 
for elements in the radioactive chain ^^Ni — > ^®Co — > ^^Fe, 
for which the intervening decay is taken into account. Con- 
sequently, only one new abundance shell is introduced, with 
boundaries given by the photospheric velocities of the pre- 
vious and the current model. The abundances in the new 
shell are set so as to achieve a best fit to the new spec- 
trum, while the previously determined abundances in the 
outer shells are retained. This procedure is carried out un- 
til all photospheric-phase spectra are modelled, yielding an 
abundance profile of the object. 

The main weaknesses of the code are the assumptions 
that the photosphere is sharp and that the radiation at the 
photosphere has the form of a black body. However, these 
assumptions are reasonable at phases up to about maximum 
light, since most 7-ray deposition occurs below the assumed 
photosphere, as discussed in Section 8. 



4 SPECTRAL ANALYSIS 



Table 1. Spectroscopic observations of SN 2002bo used for model 
calculations 



From the database presented in lBenetti et al.l i2004l) we have 
selected for modelling the spectra liste d in Table U The two 
infrared spectra already discussed in iBenetti et al.l (|20o3) 
are not included in this work. In this section the analysis of 
each spectrum is presented in chronological order. 

Global parameters are the distance, /x = 31.67, and 
the reddening. Although iBenetti et al.l i2004) favour a value 
E(B ~ V)obs ~ 0.43 ± 0.10 based on a comparison to other 
SNe la, a parameter study carried out with synthetic spec- 
tra and also presented in that paper indicates a lower value. 
The highest value that is acceptable from that study is 
E{B — V)m,od = 0.38. This value lies well within the er- 
ror range of E{B — V)obs, and is used for the calculations 
in this paper. The epochs of the spectra were derived from 



Date 


M.J.D. 


Epoch* 


Range 


Tel. 


Res. 








(days) 


(A) 




(A) 


10/03/02 


52343, 


,06 


-12.9 


3600-7700 


A1.82 


25 


10/03/02 


52343, 


,99 


-12.0 


3600-7700 


A1.82 


25 


11/03/02 


52344, 


.99 


-11.0 


3400-7700 


A1.82 


25 


13/03/02 


52346, 


.91 


-9.1 


3400-7700 


A1.82 


25 


15/03/02 


52348, 


,04 


-8.0 


3400-9050 


NOT 


14 


16/03/02 


52349, 


,93 


-6.1 


3400-9050 


NOT 


14 


18/03/02 


52351, 


.85 


-4.1 


3200-7550 


WHT 


2 


19/03/02 


52352, 


.05 


-3.9 


3400-9050 


NOT 


14 


19/03/02 


52352, 


,94 


-3.1 


3400-7700 


A1.82 


25 


20/03/02 


52353, 


,90 


-2.1 


3400-9050 


NOT 


22 


21/03/02 


52354, 


.96 


-1.0 


3400-9050 


NOT 


22 


23/03/02 


52356, 


.08 


-fO.l 


3400-10350 


A1.82 


25 


28/03/03 


52361, 


,94 


+5.9 


3100-8800 


WHT 


12 


11/12/02 


52619, 


,34 


+263 


3350-7400 


E3.6 


14 


26/03/03 


52724, 


.10 


+368 


3600-8600 


VLT 


11 



* - relative to the estimated epoch of B maximum (M JD=52356.0) 

A1.82 = Asiagol.82m telescope + AFOSC 

NOT = Nordic Optical Telescope + ALFOSC 

WHT = William Herschel Telescope + ISIS 

E3.6 = ESO 3.6m Telescope + EFOSC2 

VLT = ESO-VLT-Ul Telescope + FORSl 



Table 2. Input parameters and calculated converged temperature 
of photospheric models. 



Epoch* 


Epoch** 


Radius 


Bol. Lum. 


Temp. 






Vph 


logioL 


Tbb 


(d) 


(d) 


( km s^-"-) 


(ergs-i) 


(K) 


-12.9 


5.1 


15,800 


42.04 


9430 


-12.0 


6.0 


15,500 


42.31 


9950 


-11.0 


7.0 


15,100 


42.54 


10,930 


-9.1 


8.9 


13,900 


42.77 


11,850 


-8.0 


10.1 


12,900 


42.84 


12,330 


-6.1 


11.9 


11,450 


42.97 


13,070 


-4.1 


13.9 


10,400 


43.04 


13,310 


-3.9 


14.1 


10,200 


43.05 


13,480 


-3.1 


15.0 


9900 


43.08 


13,420 


-2.1 


15.9 


9200 


43.09 


13,820 


-1.0 


17.0 


8600 


43.10 


13,940 


+0.1 


18.1 


8100 


43.09 


13,750 


+5.9 


24.0 


7600 


43.08 


10,940 



relative to the estimated epoch of B-maximum (MJD=52356.0) 
* days after the explosion. 



a bolometric risetime of 18 d. Table|5|lists the various input 
parameters for the photospheric-epoch models, and the val- 
ues of the converged temperature at the photosphere, Tbb- 



4.1 Day -12.9 

The first spectrum in the sequence dates 12.9 days before 
maximum light, corresponding to texp ~ 5.1 d. This is one 
of the earliest SN la spectra ever taken. Fig.Qshows the ob- 
served spectrum and the best fit model. The SN is still faint, 
requiring logm L = 42.04 (ergs"'^), but the photospheric ve- 
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Figure 1. Observed spectrum and corresponding model of SN 2002bo at day —12.9 using E{B — V) = 0.38, = 31.67, teip = 5.1 d for 
the model calculation. Botii, observed spectrum and model are in observer frame. 



locity is high, Vph = 15, 800 km s ^, resulting in a low effec- 
tive temperature (Te// = 7540 K). 

Several deep absorptions dominate the spectrum and 
are labelled in Fig. Q Most of these are present through- 
out the photospheric phase. The flat overall shape of the 
spectrum a nd the ratio of the two Si ii features near 5700 
and 6OOOA iNueent et al.llT995^ are indicative of a low tem- 
perature. The feature near 7500 A (a blend of Oi 7772 A, 
Sin 7849 A and Mgll 7877, 7896 A), and the Call IR triplet 
(8498, 8542, 8662 A) near 8000 A are not blended in the syn- 
thetic spectrum, but they may be in the data. 

The near-photospheric composition includes a large 
fraction of O (30% by mass), but also remarkably high abun- 
dances of IME (including 30% Mg, 30% Si and 6% S), some 
stable Fe (1.5%), and 1.13% of '^^Ni and decay products. It 
also includes 0.02% of Ti and Cr. The high Mg abundance 
was necessary to reproduce the feature at ~ 4200 A. Ti and 
Cr are necessary to block the near-UV flux and transfer it 
to redder wavelengths. In order to reproduce the blue ex- 
tension of the Si II 6355A line an additional abundance shell 
was introduced above 22,700 km s^^. This outer shell con- 
tains mostly Si (55% by mass) but also significant amounts 
of Mg (30%), S (6%) and Ca (5.5%), so that the abun- 
dance of O is quite small. The high Ca abundance is re- 
quired to reproduce the Ca 11 lines observed at later epochs. 
SN 2002bo shows no sign of C, even at the highest velocities. 
From the present an alysis, and from that of the IR spectra 
l|Benetti et alJl2004ft . we can set an upper limit for the C 
abundance of 3% restricted to velocities > 25, 000 kms~^. C 
is not seen at later epochs either, and so we can conclude 



that essentially all C (but not O) was burned during the 
explosion. 

The main shortcoming of the model is the overestimate 
of the flux beyond about 6200 A. Since the epoch of this 
spectrum is so early that the assumption of a sharp photo- 
sphere should be correct, it is possible that data calibration 
may be the reason for the mismatch of the flux. Spectra of 
other SNe la at this early phase do not suffer from this prob- 
lem. Another possibility may be the adopted high extinction. 
However, observations suggest E{B — V)tot > 0.33, while a 
value E{B — V)tot ^ 0.25 would be needed to improve the 
situation signiflcantly. However, such a low value would lead 
to an improved model only at this particular epoch. 



4.2 Day -12.0 

The next spectrum was taken one night after the first one at 
an epoch t^xp = 6.0 d. The best-flt model is shown in Fig.|2Ji.- 
The photosphere moved inwards in velocity space and is now 
located at Vph = 15,500kms^^. The luminosity is logio-L — 
42.31 (ergs~^). Although the evolution of the spectrum is 
quite small, it did become bluer: the peaks in the U- and B- 
bands are higher with respect to those in the V- and /-bands. 
The temperature changed only marginally, as indicated by 
the Sill line ratio, whereas the model temperature increased 
by ~ 500 K, to 8150 K. This leads to a smaller Sill line ratio. 

Almost every synthetic line lacks some blue-wing ab- 
sorption. This is best seen in the Sill and Fell lines. Al- 
though the Si abundance in the outer regions is enhanced 
(see Sect. 14. ll and it is very high near the photosphere (30% 



Abundance Stratification in Type la Supernovae: I. The case of SN 2002bo 5 



3 - 



2.5 



2 - 



1 - 



0.5 - 



l"siN 
4128, 4131A 




I ' ' ' ' I ' ' ' ' I 
Day -12.0 — 
Model — 



a) 




I ' ' ' ' I ' ' ' ' I 
Day -11.0 — 
Model — 



V V .vv 



b) 





0.9 




0.8 




0.7 


•< 






0.6 


'e 










0.5 


'w 




g> 




CD 


0.4 


o 






0.3 












0.2 




0.1 









1.6 




1.4 


<< 


1.2 


OJ 




'e 




o 


1 


'cn 










0.8 






O 






0.6 


<< 




u. 






0.4 




0.2 








3000 4000 5000 6000 7000 8000 9000 

Rest Wavelength A, [ A ] 
I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
Felll Day -9.1 



4396, 4420, 4431 A 




Model -- 



c) 



I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I ■ ■ ■ ■ I 



3000 4000 5000 8000 7000 8000 9000 
Rest Wavelength A, [ A ] 
I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

Day -8.0 - 

Model 




1 - 



3000 4000 5000 6000 7000 8000 9000 
Rest Wavelength >. [A] 
I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

Day -6.1 ■ 

Model _ 




3000 4000 5000 6000 7000 8000 9000 
Rest Wavelength X [k] 
j_ I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I J. 
Ni II, Till, Day -4.1 

Co II, Co lll.i l^'lodel 
CrII, CrI.., ,1 , 

Si III 5740A 



i» 1 



3000 4000 5000 6000 7000 8000 9000 
Rest Wavelength A, [ A ] 



3000 4000 5000 6000 7000 8000 9000 
Rest Wavelength X [A] 



Figure 2. Observed spectra and corresponding models of SN 2002bo. Input values for the models are; E{B — V) = 0.38, /U = 31.67, a) 



: 6.0 d, b) te 



7.0 d, c) te 



: 8.9 d, d) te 



10.1 d, e) te 



11.9 d, f) te 



13.9 d. All models and observed spectra are in 
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by mass), this is still not sufficient to fit the blue wings of the 
lines. The line velocity required to account for the blue ab- 
sorption is > 20, OOOkms"^. To reproduce the observations 
may require an increase of the density at those velocities, 
a situ ation similar to that of most other SNe la with early 
data iMazzali et al.ll2005allbl) . 

There is still a discrepancy in the red flux between the 
model and the observations, although perhaps somewhat 
smaller than in the previous model. 

The abundances are similar to the previous epoch, but 
there is a tendency for O to decrease and for most IME 
(except Si and S) to increase near the photosphere. 

4.3 Day -11.0 

The next spectrum in the series, shown in Fig. |5j3, was ob- 
served one day later {t^xp = 7.0 d). The luminosity increased 
by almost a factor of 2 with respect to the previous epoch 
to logioL — 42.54 (ergs~^). The photosphere receded to 
Vph = 15, 100kms"\ leading to Te// = 8750 K. The ionisa- 
tion is consequently increased. The near-photospheric abun- 
dances of IME and Fe-group elements are further increased 
in this model. Although Si dominates near the photosphere, 
the synthetic Si ll lines are somewhat weak. The problem is 
that from Vph up to ~ 20, 000 km s~^ more than 99.9% of 
Si is doubly ionised. Thus, if the Si abundance is further 
increased, Sim lines appear that are not seen in the ob- 
servations. On the other hand, the strongest Sin lines are 
saturated near the photosphere and do not become stronger, 
while outer regions do not contribute since the density is too 
low. This problem will become more serious at later epochs. 

4.4 Day -9.1 

Figure 121; shows the spectrum at day —9.1 (t^xp ~ 8.9 d). 
The photosphere was located at 13,900kms~^, and the lu- 
minosity increased signiflcantly, to logioL = 42.77 (ergs~^), 
so that Tef f — 9230 K. The Si ll line ratio is well reproduced, 
indicating that the temperature is correctly evaluated. The 
weak Sill line at 4130 A is also well reproduced. 

The feature near 4200 A, which was dominated by Mg ll 
at earlier epochs, now has a signiflcant contribution from 
Fe III 4420 A. Fe is now doubly ionised over almost the entire 
envelope, which may be an overestimate since the broad Fell 
absorption near 4800 A is too weak in the model. 

The synthetic Ca ll H&K and Si ll 6355 A lines lack high 
velocity absorption. Since the abundances at higher veloci- 
ties are constrained by the previous epochs, a possible solu- 
tion is an increased density in the outermost regions of the 
ejecta, as may occur if circumstellar interaction takes place. 

The abundances continue to evolve as before, with O de- 
creasing and Fe-group elements increasing. Among the IME, 
Ca and S increase, but Mg decreases, recovering more nor- 
mal values after the large increase at the highest velocities. 

4.5 Day -8.0 

At day —8.0 {texp = 10.1 d, see Fig. |5|1) the photosphere 
reached a velocity of 12, 900 km s~^. The luminosity in- 
creased to logioi = 42.84 (erg s~^), and the effective tem- 
perature is very similar to the previous epoch (Te// — 



9390 K). The Sim line near 4400 A, which was visible in 
the models from day —11, now begins to appear in the ob- 
servations. This spectrum is of particular interest because 
it is the flrst one in the sequence that covers the Call IR 
triplet. The synthetic line appears overall somewhat weak, 
but the lack of blue absorption is particularly evident. A 
similar behaviour can be seen in the deep absorption near 
7500 A: the red part of this feature may be a little overesti- 
mated, but the blue wing is missing, probably owing to the 
weakness of the high velocity Si absorption. Other Si ll lines 
confirm this effect. 

The trend of the abundances is as at the previous epoch, 
but now both Si and Mg decrease with O. Stable Fe is not 
necessary at the velocities probed by this spectrum. 

4.6 Day -6.1 

The spectrum on day —6.1 (ieip ~ 11.9 d) is shown in 
Fig- lit- A- luminosity of logioL = 42.97 (ergs~^), and Vph = 
ll,450kms-i result in T^ff = 9860 K. This is 470 K higher 
than in the previous epoch, leading to a smaller Sill line 
ratio and to a stronger Sim line near 4400 A. The Sll ab- 
sorption near 5300 A is weaker in the model than in the ob- 
servations. This is a known problem from previous analyses 
of SNe la, possibly caused by uncertain ^/-values. 

The abundance of O and Si is unchanged, those of ^®Ni, 
S and Ca increase, compensated by the decrease of Mg. 

4.7 Day -4.1 

The spectrum at day -4.1 (Fig.|5f) has t^xp = 13.9 d. The 
model was computed with logioL ~ 43.04 (erg s~^), Vph ~ 
10,400kms"\ and it has T^ff = 10, 040 K. 

The model matches the observation very well in the 
blue, showing the ability of the MC code to treat line blan- 
keting. The absorption at the blue edge of the spectrum 
is a blend of many Niii, Coll, Com, Tin, Cm and Criii 
lines, which cannot be distinguished unambiguously. Differ- 
ences between the model and the observations are due to 
rather strong lines of Sim (~ 4400 A and 5500 A), suggest- 
ing that the temperature may be somewhat too high, and to 
missing high-velocity absorption {v > 18, OOOkms"^), most 
prominently in Si ll 6355 A and in the Mg ll-dominated fea- 
ture near 7400 A. The problem with Sll was discussed in 
Section 

Interestingly, the red flux discrepancy is much smaller 
in this WHT spectrum than in all previous ones, which were 
from Asiago and NOT. Maybe the assumption of a photo- 
spheric black body is more realistic here, but the observed 
spectrum is also signiflcantly different from all previous ones. 
Therefore, questions about the calibration remain. 

The abundance of O and Mg is further decreased, but 
other IME remain constant, having reached the maximum 
of their abundances. ^'^Ni increases further, and is the dom- 
inant element near the photosphere. 

4.8 Day -3.9 

The spectrum at day —3.9 is shown in Fig.|^. The model has 
texp = 14.1 d, Vph = 10,200kms-\ logioL = 43.05 (ergs-^) 
and Teff = 10, 080 K. 
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Figure 3. Observed spectra and corresponding models of SN 2002bo. Input values for the models are; E{B — V) = 0.38, jU = 31.67, a) 
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Although this spectrum was observed only 4.8 h after 
the one on day —4.1, it is significantly bluer. Spectral evo- 
lution is unlikely to be so rapid, highlighting the problems 
of calibration. The two spectra were observed with differ- 
ent telescopes. The exact colour evolution A(_B — V) is not 
known because no photometry is available for the epoch of 
the day —3.9 WHT spectrum, and only a _B-band measure- 
ment exists at the epoch of the NOT spectrum. The near 
photospheric composition is insignificantly different from the 
previous epoch. The model fails to reproduce the bright blue 
peaks. This may be due to the black body distribution plac- 
ing too much flux in the red (Sect. 14.6^ . but relative flux 
calibration may again be questionable. 



4.9 Day -3.1 

The spectrum at day —3.1 (teaip = 15.0 d) is shown in 
Fig. Ot. The luminosity logioL = 43.08 (erg s"^) mcreases 
more slowly compared to earlier epochs, as the light curve 
approaches its peak. The velocity decrease is also slower 
{vph = 9900 kms~^). Consequently, the temperature at the 
photosphere is the same as in the previous epoch (Te// = 
10, 090 K). The small Sill line ratio at this epoch is repro- 
duced very well by the model. 

Because of the small step in Vph, the abundances are 
practically unchanged near the photosphere, except for a 
small increase of ^®Ni and a decrease of O, Si and Ca. A 
high Si abundance is still necessary to fit the Sill 6355 A 
absorption. 



4.10 Day -2.1 

The spectrum at day —2.1, shown in Fig. |2J;, was taken 
at an epoch of texp = 15.9 d. The luminosity (logioL = 
43.09 (erg s^^)) increased only slightly with respect to 
the previous day. The photosphere receded to v^h = 
9200 km s~^, leading to a somewhat higher temperature, 
Tejf = 10, 210 K. The relative depth of the Mgll line near 
7500 A is now correctly reproduced, but both this line and 
Call IR lack blue absorption. 

The O abundance is further decreased, compensated by 
an mcrease of ^^Ni. Some stable Fe is again required, and 
all IME decrease. 



4.11 Day -1.0 

The last spectrum before maximum light has a fiducial epoch 
of t^a:p = 17.0 d (Fig.OJi). The photosph eric velocity Vph re- 
ceded to 8600 kms~^. The bolometric luminosity reached its 
maximum, logioL = 43.10 (ergs~^), at this epoch, one day 
before B maximum. The temperature also reached its max- 
imum, at a value Te// ~ 10, 250 K. The day —2.1 and the 
day —1.0 spectra look very similar, but the latter is slightly 
bluer. Unfortunately, only _B-band photometry is available 
to calibrate the spectra. A blend of Com, Felll and some 
weaker S ill lines produces the deep absorption feature near 
4300 A. The abundances of these elements may be slightly 
too high but cannot be reduced without a negative influence 
on the overall shape of the synthetic spectrum. Aundances 
have the same trend as in the previous epoch. 



4.12 Day +0.1 

The maximum light spectrum is shown in Fig.|^. An epoch 
of 18.1 days is assumed. The photospheric velocity decreased 
to 8100 kms~^, and the luminosity already begun to decline 
at logio-L = 43.09 (ergs"^). The expansion and the declining 
luminosity reduce the temperature to T^ff = 10, 200 K. 

It is again interesting that this spectrum, the only one 
extending beyond 9000 A, is also the one whose red flux is 
best reproduced by the model, even though the photospheric 
approximation is certainly less valid at maximum than at 
earlier epochs. This looks like rather convincing evidence 
that most spectra that do not extend sufficiently to the red 
are affected by significant red-end calibration problems. 

The depth of the Ca ll IR triplet relative to the contin- 
uum is reproduced very well, and only the blue wing of the 
line is missing. Si ll lines are also accurately reproduced, but 
there are unwanted Si ill lines at various places in the spec- 
trum (e.g. 4400, 5600A). The observed spectrum exhibits a 
much stronger peak ~ 3500 A than the previous one. The 
somewhat erratic relative behaviour of the blue and the red 
part of the spectrum again suggests that there may be in- 
consistencies in the data calibration. 

Abundances change smoothly, following the trend of de- 
creasing O and IME, increasing ^'^Ni and stable Fe. 

4.13 Day +5.9 

The last spectrum is from day +5.9, corresponding to t^xp ~ 
24.0 d. The observed spectrum consists of two distinct parts, 
as shown in Fig. I^T, making it difficult to tell whether the 
high emission near 5500 A is real and whether the relative 
flux level of the two parts is correct. The blue part ends just 
to the red of the S ll feature, which appears to be blended. 

Although the epoch is several days after maximum, 
the luminosity has decreased only slightly (logio-L = 
43.08 (ergs-^)). This may be a calibration problem. The 
model velocity decreased to Vph = 7600 km s~^, so the tem- 
perature dropped significantly {T^ff = 9100 K) as a conse- 
quence of the expansion. The synthetic spectrum reproduces 
the overall flux distribution and almost all of the lines. 

This spectrum explores sufficiently deep layers that the 
derived abundance distribution overlaps with that from the 
models of the nebular phase (Section [KJ. However, the as- 
sumption of a sharp photosphere is weakest at this epoch, 
making the derived abundances more uncertain. Neverthe- 
less, the fact that the abundances continue to change in the 
same way as at earlier epochs is comforting. 



5 NEBULAR PHASE 

The photospheric epoch coverage of SN 2002bo ends at day 
+6, when Vph ~ 7600 km s~^. It is therefore not possible to 
extend the analysis to lower velocities using the MC code. 
In any case, the code becomes less and less applicable at ad- 
vanced post-maximum epochs, as the photospheric approx- 
imation becomes rapidly inadequate. There is, however, an 
alternative approach that makes it possible to investigate 
the properties of the inner part of the ejecta. That is mod- 
elling the spectra in the so-called nebular phase, when the 
ejecta are completely transparent. Models applicable to this 




Figure 4. Observed spectra and corresponding models of SN 2002bo in tiie nebular phase. Input values are: E(B — V) = 0.38, fi = 31.67, 
a) texp = 282 d, b) texp = 386 d. 



epoch ha ve been developed, an d applied to various types of 
SNe fe.g. lMazzali et alji2001bh . These models follow the en- 
ergy deposition and heating by the 7-rays and the positrons 
from the radioactive decay of ^®Ni and ^^Co, and the ensuing 
cooling by line emission. 

We extended our previous one-zone model also to treat 
abundance stratification. Th e 7- ray deposition is treated in 
a Montecarlo approach (see ICap pcUaro ct al. 1997) , while 
nebular emission is computed in NLTE. This code can be 
used to reproduce accurately the emission line profiles, and 
hence to derive the abundance distribution in the ejecta. 

Two nebular spectra are available, taken on Dec. 11th, 
2002 and March 26th, 2003, i.e. 282 and 386 days after the 
estimated explosion epoch, respectively (see Tab.0. They 
were reduced with standard IRAF^ routines. Extractions 
were weighted by the variance based on the data values and 
a Poisson/CCD model, using the gain and read noise param- 
eters. The background to either side of the SN signal was fit- 
ted with a low-order polynomial and then subtracted. Flux 
calibration was performed using spectrophotometric stan- 
dard stars. 

We modelled the spectra adopting the W7 density pro- 
file, fixing the abundances in the outer part of the ejecta 
according to the results of the photospheric epoch study 
above. Thus we only had the freedom to modify abundances 
below 7600kms"\ 

The model for the day -1-263 spectrum is shown in 
Fig. UK- It reproduces the observed spectrum very well if 
the abundance of ^^Ni at velocities below 7600 kms~^ is in- 
creased above the values estimated from the photospheric 
epoch study at higher velocities. However, the ^®Ni abun- 

^ The Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by the National Optical Astronomy Observatory, which 
is operated by AURA, Inc., under a cooperative agreement with 
the National Science Foundation. 



dance reaches at most 0.62 by mass in regions between 4000 
and 7000 km s~^. This confirms that significant mixing out 
of ^®Ni occurred, as derived from the early-time spectra. Fur- 
ther inwards, the ^"^Ni abundance decreases, as expected. If 
the innermost ejecta contain mostly ^®Ni, the Fe-line emis- 
sion becomes too strong. However, line profiles show that 
there is Fe emission at the lowest velocities. This can be re- 
produced if the innermost layers contain significant amounts 
of stable Fe (~ 0.30 Mq at v < 5000 kms"^). Stable Fe 
must be excited by (mostly) 7-rays and (almost negligibly) 
positrons emitted at higher velocity layers and difi^using in- 
wards. Stable Fe should be mostly ^*Fe and ^*Fe. These 
isotopes are the result of burning to NSE, but at slightly 
different densities. Spectroscopically, we are unable to dis- 
tinguish between them. A total ^^Ni mass of 0.50 Mq was 
used, which makes SN 2002bo an average SN la. The '^^^Ni 
abundance is low. A total ~ 0.80 Mq was burned to NSE. 
The abundance distribution is plotted in Fig. 5. 

Synthetic spectra for the later epoch (day -1-368) com- 
puted using the same abundance distribution nicely confirm 
these results (Fig.QJj). 

Significant changes of the line profiles or width over time 
are not visible in the data. 



6 ABUNDANCE STRATIFICATION 

Modelling of all spectra, both of the photospheric and the 
nebular phase, delivers the abundance distribution of the en- 
tire ejecta. The relative abundance by mass of the elements 
that could be detected in the spectra is shown in Fig. as 
a function of radius in velocity space. 

The various steps above 7600 km s~^ reflect the radii 
of the early-time spectra. Therefore, the distance between 
two abundance shells depends on the time interval between 
the two spectra. An additional shell above 22, 700 km s~^ 
was inserted in order to account for high velocity material 
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Table 3. Total mass of elements spectroscopically confirmed in 
the ejecta of SN 2002bo and nucleosynthesis products of different 
SN la explosion models 



Ejected (Synthesized) Mass (Mq) 



Species 


SN 2002bo 




WDDf'' 


b30.3d.768'= 


C 


< 0.002 


4.83E-02 


5.42E-03 


2.78E-01 


O 


O.flO 


f.43E-0f 


8.82E-02 


3.39E-0f 


Na 


O.OOf 


6.32E-05 


8.77E-05 


8.65E-04 


Mg 


0.080 


8.58E-03 


7.69E-03 


8.22E-03 


Si 


0.220 


f.57E-0f 


2.74E-01 


5.53E-02 


S 


0.067 


8.70E-02 


1.63E-01 


2.74E-02 


Ca 


020 


f.l9E-02 


3.fOE-02 


3.61E-03 


Ti 


0.003 


3.43E-04 


f.f3E-03 


8.98E-05 


Cr 


0.003 


8.48E-03 


2.05E-02 


3.f9E-03 


Fe* 


0.360 


f.63E-0f 


1.08E-01 


l.f3E-01 


56^;** 


0.520 


5.86E-0f 


5.64E-0f 


4.f8E-0f 


Ni*** 


n.a. 


1.26E-0f 


3.82E-02 


f.06E-0f 



* Stable isotopes except for ^^Fe from ^^Co decay 
**Isotopes from the ^^Ni — > ^^Co — > ^^Fe decay chain 
" ""Stable isotopes ^^N i, ^°m, "^^Ni, 62Ni, and ^"Ni 
° Nomoto et alJ il984h 
Iwamoto et alj il99'gi) 
'' iTravaglio et alj \2Q(m 



(Sect. 4.1). Below 7600 kms~^, where the abundances are de- 
rived from the nebular spectra, the radii of the shells are 
given by the density shells of the underlying W7 model. 
Since our models cover the entire SN ejecta, the total mass 
of each element can be calculated. These values are listed in 
Table 13 Different isotopes cannot be distinguished spectro- 
scopically, so here the abundance of an element implies the 
sum of all isotopes, except for Fe. 

Since the progenitor WD is supposed to be composed 
of C and O, the abundances of these two elements in the SN 
ejecta indirectly indicate the extent of burning. Interestingly, 
no sign of C can be seen in the spectra of SN 2002bo. The 
most prominent C lines in the optical and NIR wavelength 
bands are C ii 6579 A a nd C ii 7231 A Neith er hne can be 
detected, at any epoch. iBenetti et al.l i2004l) suggested an 
upper limit of the C abundance of 3% at v > 25, OOOkms"^ 
based on models of the optical spectra. This upper limit 
corresponds to a total mass of < 2 X 10"^ Mq. Since no C 
lines are detected in the spectra, the C abundance in this 
analysis is set to zero. 

The other progenitor element, O, dominates the ejecta 
between 15, 500 kms~^ and 22, 700 kms~^. In the outermost 
layers the abundance of Si is higher than that of O. In deeper 
layers the O abundance decreases, and no O can be found 
below 6000 kms"^ Altogether we detect ~ 0.11 M© of O in 
SN 2002bo. 

The group of IME includes Na, Mg, Si, S and Ca. Ex- 
cept for Na, they all show rather high abundances at high 
velocities. Na itself is located above the Fe core and out to 
~ 8500 km s~^, with a low abundance. 

Mg is enhanced in SN 2002bo. A total mass of ~ 
0.08 Mq is derived from the modelling. The bulk of Mg is at 
V > 12, 000 km s~^, with relative mass abundances between 
20% and 30%. This is necessary to fit the deep 4481 A ab- 
sorption in the early phase, as discussed in Section |1] Below 



~ 15, 000 km s~^ the Mg abundance decreases steadily, and 
no Mg is detected below 7300 km s~^. 

The Si abundance has a similar behaviour as Mg. How- 
ever, the increase of the abundance at high velocity is more 
pronounced. The outermost shell (« > 22, 700kms~^) con- 
tains 52% of Si in order to account for the blue wing ab- 
sorption of the various Si ll features. In the region between 
11, 450 km s"^ and 22, 700 km s"^ Si is as abundant as O 
(~ 30%). This is similar to the prediction of a model like W7, 
although somewhat on the high side. Deeper in the ejecta 
the Si abundance is ~ 15% down to ~ 6000 km s~^, and then 
it drops rapidly until it disappears below 4000 kms^^. A to- 
tal mass of ~ 0.22 M0 is derived, which is more than in W7 
but less than in the DD models. It must be emphasised that 
the Si abundance in the deepest layers can only be derived 
indirectly, since the only nebular lines of Si are in the IR, 
which is not available for SN 2002bo. 

S behaves like Si, but with a lower abundance, and 
does not exhibit a rise in the outermost shell. It starts 
at 6% outside, increasing to ~ 13% between 12,900 and 
10,200 km s~^, and decreasing to ~ 6% before disappearing 
below ~ 6000kms"\ The total mass of S is ~ 0.067 Mq. 

Ca completes the set of IME that can be detected in 
SN 2002bo. This element is the best to study high velocity 
components, since its lines are the strongest. As discussed 
in Sect. 2] the Call IR triplet lacks significant absorption 
in its blue wing. We tried to account for this using up to 
20% Ca at intermediate and especially at high velocities. 
However, even with 100% Ca in the outermost layers it is 
impossible to fit the blue part of the line in the day -8 
spectrum. On the other hand, the strong Call 7291, 7324 A 
nebular emission can be reproduced with a very small Ca 
abundance (~ 10~^ by mass). Therefore, a second iteration 
of the modelling of the early time spectra was performed 
using a significantly lower value for the Ca abundance. A 
satisfactory result was obtained for an abundance of 5.5% 
above 22, 700kms~^, decreasing to 0.3% in the next deeper 
shell and increasing again to 1% - 5.5% between 15, 800 
and 7500 km s~^. Below this velocity, the abundance is very 
small. A total 0.021 Mq of Ca is estimated. The problem 
with the high velocity absorptions cannot be solved by a 
simple enhancement of the abundances, and may require 
scenarios like an increased density in the outer ejecta, pos- 
sibly due to circumstellar interaction. 

For Ti and Cr, which have a similar abundance trend, 
a total mass of ~ 6 x 10~^ Mq is derived. Although no indi- 
vidual lines of these species can be detected in the spectra, 
their line-blocking effect is necessary to shift UV and blue 
photons to longer wavelengths. The distribution of these ele- 
ments peaks between 10, 000 and 15, 000 km s^^ , with values 
between 1 and 2%. 

The abundance of ^®Ni includes plus both the ^^Co and 
^^Fe that form in the decay chain. The relative abundances 
of these species change with time. ^®Ni dominates between 
about 3,000 and 10000 kms~^, but its abundance never ex- 
ceeds ~ 60%. Beyond this velocity the abundance decreases, 
dropping to < 1% above 15, 500kms~^. We estimate the to- 
tal 5«Ni mass synthesised to be ~ 0.52 Mq. 

The abundance of just Ni includes all stable Ni nuclei 
(^**Ni, ^°Ni, ''^Ni, ''^Ni, and "^Ni), which are mostly located 
at the lowest velocities. The dominant isotope should be 
^*Ni. The production of these species depends on the neu- 
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tron excess during the burning regime. Especially in the W7 
model stable Ni makes up a significant mass fraction of the 
inner ejecta. The synthetic spectra cannot differentiate be- 
tween different isotopes. 

The heading "Stable Fe" stands for the sum of all Fe iso- 
topes C^^Fe, '^'^Fe, ^''Fe, ^*Fe), except for the '^'^Fe that is pro- 
duced in the decay of ^®Co. A stable Fe core extends out to 
~ 3000 km s~^. Further out stable Fe decreases steadily, and 
it is absent in the velocity range from 9900 to 12, 900kms~^. 
This means that the ^''Fe from the decay of ^''Ni is sufficient 
to reproduce the spectra. Above this void region a small 
amount of stable Fe (~ 2%) is again detected, extending to 
the highest velocities. The presence of stable Fe at high ve- 
locities is required because the outer layers are observed so 
early that only very little *^Fe has been produced from ^®Ni. 
Altogether a total mass of 0.36 of stable Fe is measured. 

We compared the derived abundance distribution of 
SN 2002bo with various explosion models. Figure [KJj shows 
the distribution of the main elements (O, Mg, Si, stable Fe 
isot opes, and ^^Ni) versu s M(r). The 1-D deflagration model 
W7 llNomoto et alll984l) is shown in Fig.|3; for comparison. 
The radial distribution is very similar in the two plots. In 
both cases, the innermost part of the envelope is dominated 
by stable Fe isotopes. Above this is the ^^Ni-dominated re- 
gion, and IME are located further out. Despite these sim- 
ilarities there are significant differences. In SN 2002bo the 
abundance pattern is shifted to higher velocities and the re- 
gions of Fe group elements, IME, and unburned material 
are not as well separated but rather more mixed. The inner 
stable Fe region is much larger. O is mixed down. Si and 
Mg extend to both lower and higher velocities than in W7 
while ^^Ni and Fe are mixed to the outer regions. Finally, 
the region above ~ 1.3 Mq, which is completely unburned 
in W7, in SN 2002bo shows no trace of C and both IME 
and Fe-group elements are mixed outwards. In practice, in 
SN 2002bo we detect more stable Fe but no stable Ni (mostly 
'5**NiinW7). 

Since delayed detonation models burn high velocity re- 
gions to IME, they could possibly compare better to the 
abundance distribution of SN 2002bo. Hence we examined 
various detonation models , namely WDDl, WD D2, WDD3, 
CDDl, CDD2, and CDD3 hwamoto et alll99g^ . Indeed, the 
occurrence of IME at high velocities in SN 2002bo is repre- 
sented better by the DD models than by W7. However, O 
and the Fe-group elements show a very different behaviour, 
much closer to W7. In particular, the relative locations of 
Fe, Ni, IME, and unburned O is reproduced best by W7. 

Therefore we may conclude that in the case of 
SN 2002bo W7 is the explosion model that fits best the 
abundance distribution derived from spectral synthesis cal- 
culations. However, if mixing effects are neglected it is dif- 
ficult to reproduce the observed abundance distribution ex- 
actly. On the other hand, the analysis of a single object is 
not sufficient to constrain different explosion models. Abun- 
dance stratification analysis of several objects is necessary, 
and it is under way. 



7 BOLOMETRIC LIGHT CURVE 

In order to check the validity of our results independently, we 
used the abundance distribution obtained from the spectral 
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Figure 6. uvoir light curve of SN 2002bo. The black dots repre- 
sent the observed photometric points, the circles are the bolomet- 
ric luminosity deduced from the synthetic spectra. Three models 
of the MC light curve code are shown using directly the abun- 
dances from the spectral analysis ( ), the same abundances 

but with slightly increased ^^Ni mass ( — ), and W7 abundances 
with 0.5 M0 of 56 Ni ( ). 

modelling and computed a synthetic bolometric light curve, 
assuming that the de nsity is given by W7 and using a grey 
MC lig ht curve model ijCaopellaro et al.ll997l : lMazzali et alJ 
[20^^. In Fig.|S|we compare the resul ts with the bolometri c 
light curve of SN 2002bo presented in lBenetti et alJ J2004 . 

The synthetic light curve reproduces the observed one 
successfully, and it performs exceedingly well in the early, 
rising branch. This is particularly exciting. In fact, the early 
light curve rise depends mostly on the mixing out of ^"^Ni. 
A synthetic light curve computed with the unmixed W7 
abundances and scaled to a ^^Ni mass of 0.5 Mq shows a 
significantly later rise, which is in disagreement with the 
observations (Fig.jSJ. In that case, ^^Ni extends only up to 
~ 20, 000 km s~^ , so photons are trapped longer in the ejecta 
and the rise of the light curve is delayed. Evidently, abun- 
dance stratification is accurately determined by the spectral 
modelling, especially in the outer layers. 

The model is less successful near the peak of the light 
curve. A slight modification of the abundances, with ^^Ni 
increased to about 0.56 M0 and stable Fe decreased accord- 
ingly, yields a better fit of this phase. This is a small change, 
which does not affect the essence of our results. 



8 THE POSITION OF THE PHOTOSPHERE 
AND THE RELIABILITY OF THE RESULTS 

We mentioned that the velocities of IME lines in SN 2002bo 
are peculiarly high. Therefore, it is interesting to see how 
the velocities of the observed lines compare to the model cal- 
culations. Especially for the weaker lines, the photospheric 
radius should be a good estimate for the expansion velocities 
of the model line features. The evolution of the photospheric 
velocity is compared in Fig. |7|to the expansion velocities of 
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Figure 7. Evolution of the expansion velocities deduced from the 
minima of the Si II 6355 A, S II 5640 A, Call H&K absorptions and 
the photospheric radius of the synthetic spectra, and the position 
of the photosphere in the light curve code. 



Sill 6355 A, S ll 5640 A, and Can H&K, as measured from 
their minima llBenetti et al.ll2004^ . 

The photosphere of SN 2002bo is at 15,800 km s"^ on 
day —12.9, and recedes only by 700kms~^ until day —11.0. 
This slow decline phase accounts for trying to reproduce 
the high- velocity absorption features. After day —11.0, the 
decUne is steeper (590kms~^ d~^), and then it slows down 
again after maximum. The three lines exhibit a slightly dif- 
ferent behaviour than the photosphere. Their velocities de- 
crease more rapidly until day ~ —8 then they do afterwards. 
After maximum, the gradient is very similar to the photo- 
spheric one. 

The S II line has the lowest velocity, as may be expected 
since this is the weakest line that can be individually mea- 
sured. After day —8.0 it forms well above the Montecarlo 
photosphere, but it matches the light curve photosphere up 
to abouty maximum. 

The strong Si ll line declines at a slower rate, especially 
between day —6.1 and maximum. It is always significantly 
faster than the photospheric velocity, as can be expected 
since this strong line is formed well above the photosphere. 

The Ca ll H&K line is very strong and it has the highest 
velocities. It evolves like the Sill line, at higher velocities. 
However, the very high velocities measured in the earliest 
spectra (~ 25, 000 km s"^) in both Call H&K and Call IR 
are probably not photospheric, but indicate the presence of 
a high- velocity component (Mazzali ct al. 2005a b). 

The calculation of the light curve offers the opportunity 
to investigate the validity of the assumption of a sharp pho- 
tosphere that is made in the MC code. The light curve code 
in fact computes 7-ray deposition as a function of depth 
and time, and estimates the position of a gray photosphere 
based on a simple opacity prescription. The photospheric 
velocity thus computed, also shown in Fig. 6, follows closely 
the position determined by the MC spectrum synthesis code, 
and, even more closely, that traced by the velocity of the S 11 



line, which is sufficiently weak to be a good tracer of Vph- We 
can now look at the fraction of the energy deposited above 
the photosphere as a function of time. This fraction reaches 
10% between 1 and 2 days before maximum, it is ~ 13% 
at maximum and ~ 25% 6 days after maximum. These val- 
ues indicate that up to about maximum the approximation 
of a sharp photosphere is acceptable, while in the spectrum 
at day +5.9 it is rather poor, and non-thermal excitation 
and ionisation effects that are not considered in the code 
may be relevant. Another test is the fraction of ^^Ni located 
above the photosphere. This is only ~ 10% at maximum and 
~ 20% 6 days later, confirming the above result. 

Fortunately, the depth of the zone explored between 
maximum light and the nebular spectra is quite small (be- 
tween 7600 and 9000 km s~^), since the photosphere recedes 
only slowly after maximum and then disappears. Neverthe- 
less, abundances in that velocity range must be regarded as 
more uncertain than values both at higher velocities, where 
the MC code is more reliable, and at lower ones, which are 
treated with the nebular code. 



9 CONCLUSIONS 

We have derived the abundance stratification in the ejecta of 
the SN la 2002bo through a time series of spectral models. 
Synthetic spectra were computed for 13 epochs during the 
photospheric phase and 2 in the nebular phase. 

The most important result is that although the abun- 
dances are not very different from those of a standard model 
such as W7, some mixing in abundance seems to have oc- 
curred. In particular, ^®Ni extends to higher velocities than 
in W7, and the IME are mixed both inwards and outwards. 
Some IME are present at the highest velocities, which may 
indicate interaction of the ejecta with a CSM or some as- 
phericity in the explosion. Balancing this, O is mixed in- 
wards. C is however not seen in the ejecta. 

One possibility to produce a situation like that observed 
in SN 2002bo is multi-dimensional effects. Bubbles of burned 
material may rise to outer regions while fingers of unburned 
material sink to deeper layers. Unlike the case of global mix- 
ing, the direction from which the observer looks at the SN 
may be important. In the case of SN 2002bo one of these par- 
tially burned bubbles reaching out to high velocities may be 
in the line of sight, delivering the spectra and abundance 
distribution we observe. Had we looked from another posi- 
tion we might possibly have seen the missing C and less or 
no high-velocity c omponents in the lines of IME (see, e.g. 
iKasen et al.ll200'3) . In this context it would be interesting to 
verify whether multi-dimensional effects can account for the 
observed diversity in SNe la beyond the brightness-decline 
rate relation. 

The investigation of the abundance stratification in a 
number of well-observed SNe la can contribute to our un- 
derstanding of all these mechanisms. 
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